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Abstract 

The scaled-particle theory equation of state for the two-dimensional 
hard-disk fluid on a curved surface is proposed and used to determine 
the saddle-splay modulus of a particle-laden fluid interface. The resulting 
contribution to saddle-splay modulus, which is caused by thermal motion 
of the adsorbed particles, is comparable in magnitude with the saddle- 
splay modulus of a simple fluid interface. 



1 Introduction 

The surface free energy density of fluid interfaces depends upon their curvature. 
This dependence affects the nucleation in Uquids [TJ [51 [3] , and has important 
role in determining the structure and dynamics of the systems with complex 
fluid interfaces, such as membranes or surfactants [1]. 

For small curvature of the interface, the dependence of the surface free energy 
/ upon the geometry of the interface is conveniently described by the Helfrich 
curvature expansion [5] 

J = a + 2k{H - H^f + RK (1) 

In this equation the geometry of the interface is characterized by mean curvature 
H — + I/R2) and Gaussian curvature K = l/(i?i_R2), ^^i and R2 being 

principal radii of curvature of the interface. The surface tension a, bending 
modulus K, spontaneous curvature Hq, and saddle-splay modulus (or Gaussian 
rigidity) R, are the material parameters of the interface. By virtue of Gauss- 
Bonnet theorem, the contribution of the last term in Eq. ([1]) into the total free 
energy of the system depends on the topology of the system. Indeed, the value 
of the saddle-splay modulus affects the processes which involve changes in the 
topology of fluid interfaces H [HI El [IHl E] . 

An interesting example of a system which can be macroscopically viewed 
as a complex fluid interface is the fluid interface laden with colloidal micro- or 
nanoparticles. To minimize total interfacial energy, particles suspended in a 



1 



bulk fluid self-assemble on the fluid interface [12] • This process, flrst observed 
by Ramsden in 1903 [13j, has recently attracted significant scientific attention 
[HI [T51 [TBI [TT] . It has also potential for a range of novel applications [IHl dSl 

On the scale large compared to the size of the adsorbed particles, a particle- 
laden interface may be viewed as continuous. If the interface is isotropic on this 
scale, the interfacial free energy can be described by Eq. ([1]), and the interface 
can be characterized by the material parameters a, k, Hq, and S. 

The present letter is devoted to the study of the saddle-splay modulus R 
of a particle-laden fluid interface at low surface concentration of the adsorbed 
particles. In this case we can represent the interface as a two-dimensional fluid 
on a curved surface. The main contribution to the interaction between particles 
at low concentration comes from the excluded volume (different particles cannot 
occupy the same space). Hence we approximate the system by a two-dimensional 
hard-disk fluid on a curved surface. 

Hard-disk fluids in curved geometry were used before to study packing of 
disks [531 [Ml US US] , ordering phase transition , topological defects [28j , and 
as a model of glass-forming liquids [M] [25l [29] . Several equations of state were 
proposed for hard-disk fluids in spherical [5D1 [ST] and hyperbolic [321 dSl 133] 
geometries. 

In the present work we shall use scaled-particle theory (SPT) [31] to derive 
the equation of state of two-dimensional hard-disk fluid on a curved surface. We 
shall then use the resulting equation of state to determine saddle-splay modulus 
R for particle-laden fluid interface at low concentration of the adsorbed particles. 



2 Saddle-splay modulus 

In accordance with Eq. ([T]), the saddle-splay modulus is given by the derivative 
of the surface free energy density with respect to Gaussian curvature, 



-_ df 



(2) 

K=0 



Using the expression for the excess free energy 

^^cx fP z -1 

N ' 

where 



dp, (3) 



(4) 

P 

is the compressibility factor, p = N/A is the number density (number of particles 
per unit area), P is pressure, /3 = l/fc^T is the inverse temperature, we may 
represent Eq. ^ in the form 



dp, (5) 

K=0 
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where the derivative is taken at constant particle density p. 

Equation ([5]) can be used to calculate saddle-splay modulus of the interface 
from the curvature dependence of the compressibility factor, which is generally 
given by the equation of state of the system. We shall use the SPT equation of 
state for a hard-disk fluid on a curved surface, which is derived in the following 
sections. 



3 SPT equation of state for hard disks 

Scaled-particle theory was originally developed by Reiss et al [34] and further 
improved afterwards [351 ISZl I3H1 132] ■ Applied to the case of hard disks 
on a 2D plane, SPT leads to a particularly simple equation of state which is 
nevertheless in good agreement with computer simulation results throughout 
most of the fluid range of densities [101 SI] • 

SPT for two-dimensional hard-disk fluids in its simplest form can be sum- 
marized as follows (see textbook [32] for more details). The reversible work 
W(i?o) is considered which is required to create a circular cavity of radius i?o 
in the fluid of hard disks of radius R. The assumption is made that for Rq > 0, 
W{Ro) is given by a polynomial in Rq 

W{Ro)^wo+wiRo + SiRo)P, i?o > 0. (6) 

The last term S{Ro)P {S{R) being the area of the disk of radius R), which is 
dominant for large cavities {Rq 3> R), follows from thermodynamics. For small 
cavities {0 < R + Rq < R), T4^(i?o) can be written in form 

WiRo) ^ -keTln [I - pS{Rq + R)] , -i? < i?o < 0. (7) 

The coefficients wq and wi are then determined by requiring the work W{R{)) 
and its derivative T4^'(i?o), given by Eqs (O and ([7]), to be continuous at i?o = 0. 
The explicit expression for the excess chemical potential of the fluid, fi"^^ = 
W{R), can be determined from Eq. and subsequently used to write the 
SPT equation of state. 

In the case of the flat surface, the area of the disk is 

S{R) = ttR^. (8) 

The corresponding values of the coefficients Wi are given by 

(3w^ = -Hl-il), /3u;i = ^, (9) 

1—7/ 

where 77 — irR^p is the hard-disk packing fraction. The chemical potential of 
the fluid, /i, is given by 

/3/z = lnAV-ln(l-r7) + ^^ + ^, (10) 

1-77 p 
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where A is de Broglie thermal wavelength. The SPT equation of state is then 
obtained from Eq. (jlOp and the thermodynamic relation 

dP dfi 

and has the form reported by Helfand et al | 40| : 

4 SPT equation of state for hard disks on a 
curved surface 

SPT equation of state for a hard-disk fluid on a curved surface can be obtained 
in the same way as in the flat case described above. The difference is that 
expression ([8|) for the area of the disk of radius R is no longer valid on a curved 
surface. For small Gaussian curvature {K <^ we shall replace it by the 

formula for the area of a geodesic disk on two-dimensional Riemannian manifold, 
obtained by Bertrand and Diguet in 1848 [13], 

5(i?) = 7ri?2(l-0 + o(0, (13) 
where we have introduced the dimensionless quantity 

Requiring work W{Ro) and its derivative W'{Ro)^ as given by Eqs ^ and ([7]), to 
be continuous at i?o = 0, we obtain the following expressions for the coefficients 

/3w;o = -ln[l-^i?V(l-e)], (15) 

and the chemical potential, 

= lnA2p-ln[l- 77(1-0] 



1-77(1-0 P 



(17) 



Equations (fT7|) and pTjl lead to the following form of the SPT equation of state 
for hard-disk fluid on a curved surface: 



^= M „./%M2 - (18) 



[1-77(1 -OF 

Figure [T] demonstrates satisfactory agreement of the compressibility factor 
Z calculated from Eq. (|18p with the Monte Carlo results for hard disks on a 
sphere reported by Giarritta et al [27] . In the case of zero Gaussian curvature 
(C = 0) Eq. (dlD coincides with Eq. 
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Figure 1: Compressibility factor Z as a function of reduced number density 
of the fluid pa^, where a = 2R is particle diameter. Circles represent Monte 
Carlo results for N — 400 hard disks on a sphere [27l , line corresponds to SPT 
equation of state (fT8|) . 



5 Saddle-splay modulus from SPT equation of 
state 

The expression for saddle-splay modulus is obtained by substituting the com- 
pressibility factor given by the equation of state, Eq. (fTS]) . into formula ([5]). The 
result is 



Note that although using the truncated series in R given by the formula ([Tc 
for the area of the large disk in the expression ^ is generally not justified, it 
is still suitable for our purpose of calculating saddle-splay modulus since we are 
interested in the limit K 0. 

The dependence of the saddle-splay modulus R upon the disk packing fraction 
7y, given by Eq. (|19p . is presented in Figure [21 The value of the saddle-splay 
modulus for particle-laden interfaces appears to be smaller than the values |k| ~ 
lOfc^T typical for lipid monolayers [IJ), but is comparable to the value |k| w 
^ksT for the surfaces of simple fluids dSj. 
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Figure 2: The dependence of the saddle-splay modulus k upon disk packing 
fraction ij, calculated using Eq. ^9]) . 

6 Conclusion 

The main message of this letter is that the thermal motion of the particles ad- 
sorbed on a fluid interface contributes to the saddle-splay modulus of the inter- 
face. This result may have implications in the structure and dynamics particle- 
laden systems that allow topological changes, for example, fusion of particles in 
Pickering emulsions, or structural reorganization in particle-stabilized foams. 

The simplest version of the scaled-particle theory allows construction of a 
rather simple equation of state for a hard-disk fluid on a curved surface. In order 
to improve the formula obtained for the saddle-splay modulus of a particle-laden 
fluid interface, it seems reasonable to attempt to construct the equation of state 
that gives more accurate dependence of the compressibility factor with respect 
to Gaussian curvature of the interface, which can be verified by using the virial 
expansion on the curved surface or the computer modelling of the system. 

The result can also be extended by taking into account the influence on 
the value of saddle-splay modulus of other contributions to the interparticle 
interaction, such as capillary, electrostatic, van der Waals etc, as well as the 
role of particles' anisotropy. The prediction of the elastic properties of the 
interfaces with large concentration of particles, in which two-dimensional solid 
structure forms, presents another interesting and more complicated problem. 
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